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The prospect of developing magnetic qubits is discussed. The first part of the article makes 
suggestions on how to achieve the coherent quantum superposition of spin states in small ferro- 
magnetic clusters, weakly uncompensated antiferromagnetic clusters, and magnetic molecules. The 
second part of the article deals with mechanisms of decoherence expected in magnetic systems. 
Main decohering effects are coming from nuclear spins and magnetic fields. They can be reduced 
by isotopic purification and superconducting shielding. In that case the time reversal symmetry of 
spin Hamiltonians makes spin-phonon coupling ineffective in destroying quantum coherence. 
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I. INTRODUCTION 

Is quantum computing physics or engineering? Some time ago I asked this question to Murray Gell-Mann who 
replied that he had a formal proof that it is physics: articles on quantum computing are published by The Physical 
Review Letters. In the spirit of his answer I will concentrate on the physics aspects of the problem, leaving out the 
engineering aspects. The purpose of this article is to discuss the prospect of using small magnetic clusters as qubits. 

Qubits based upon quantum superposition of the | t»> and | J,> states of individual electrons and nuclei have been 
discussed during all years of quantum computermaniafl. In parallel, but with little overlap, there has been intensive 
theoretical and experimental research on spin tunneling between the | | > and | |> states of molecular magnets and 
in ferromagnetic and antiferromagnetic nanoparticlescl. These are composite objects with total spin S ranging from 
S — 10 for Fe-8 and Mn-12 molecular clusters to S of a few thousand in nanoparticles. Experimental study of the 
magnetization reversal in individual nanoparticles gives evidence that very small particles are uniformly magnetized!! 
At low enough temperature they, like molecular magnets, possess a large fixed- length spin formed by the strong 
exchange interaction. 

Quantum tunneling bctsteen the | |> and | |> states of spin-10 molecular nanomagnets has been unambiguously 
established in experimental]. Early measurements of low temperature magnetic relaxationa and more recent mea- 
surements of individual nanoparticles^ also provided evidence of quantum tunneling of spin. There is experimental 
evidence of quantum coherent oscillations between the | ]> and | |> states in nanoparticlesu and molecular magnetsa. 

On one hand, individual nanoparticles and high-spin molecules, because of their large magnetic moments, must be 
easier to operate as qubits than individual electron and nuclear spins. In fact, the techniques of measuring individual 
particles of S^IO 2 — 10 3 already existB. On the other hand, as the size of the system increases, its interaction with 
the dissipative environment also increases and it is not obvious whether the decoherence in large spin systems can be 
made low enough to allow their application as qubits. To achieve this goal the decoherence rate 7 must be made small 
compared to the frequency A/ft, of coherent oscillations between the | |> and | |> states. For a large spin system it 
requires effort for two reasons. Firstly, the interaction with the environment is proportional to the size of the system. 
Secondly, log(A) scales linearly with —S so that, except for some special cases, A becomes immeasurably small for 
S>30. 

In this article we will make suggestions on how to obtain large A for a relatively large spin and how to achieve 
the condition 7 << A/ ft. It is our believe that this is doable in spin systems, in part due to their special properties 
with respect to time reversal. Spin Hamiltonians and tunneling rates will be considered in Section 2. Mechanisms of 
decoherence will be discussed in Section 3. 



II. COHERENCE 

A. Ferromagnetic clusters 

The basic Hamiltonian for a ferromagnetic cluster of spin S that exhibits quantum coherence is 

H = -AS\ + V , (1) 
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where A is a positive constant and V is an operator that does not commute with S z but is invariant under the 
transformation S z —* — S z . The first term in Eq.(l) is typically produced by a crystal field or by the shape anisotropy 
of the cluster. At V = the Hamiltonian (1) has a double degenerate ground state corresponding to two opposite 
orientations, | j> and | J,>, of S along the Z-axis. In terms of the magnetic quantum number, S z \m > = m\m >, 
these ground states are \S > and | — S >. Note that this degeneracy is independent of any geometry of the problem 
including the shape of the magnetic cluster and is solely due to the odd symmetry of S with respect to time reversal. 

Our case of interest will be an integer large S. In that case the problem is almost classical so that even a large 
perturbation Vy^O generates only a small probability of tunneling between the | |> and | |> states. The degeneracy 
of the ground state is then removed and the new ground state can be approximated by 
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It is separated from the first excited state, 
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by the energy gap A which is determined by the strength and the nature of V and is small compared to the scale 
A of the energy levels of —AS Z . The probability of finding S looking up (down) then oscillates in time according 
cos(A-t/Ti), which is quantum coherence. 

One of the possible forms of V is V = —gfJ,BH x S x induced by the external field applied along the X-axis. If the 
field is small, the tunneling splitting can be obtained by the perturbation theorjJ13 
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^From practical point of view this case is not very promising since the inevitably present weak misorientation of the 
field, resulting in H z ^0, will destroy the coherence. 

A more promising case, which to a good approximation corresponds to the Fe-8 spin-10 molecular nanomagnet, is 
V = BS X . At |B| << A the perturbation theory givesEl 
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In Fe-8 A/h is of order 10 4 s _1 . 

For large S and arbitrary B the tunneliag|SKlitting has been computed by the instanton methodeja and by mapping 
the spin problem onto a particle problenOllil 
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^From practical point of view, in ferromagnetic nanoparticles with S >> 1 the case of interest is \B\ » A. In that 
case Eq.(6) gives 



A = 167r- 1 / 2 S 3 / 2 A 3 / 4 |5| 1 / 4 exp 
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and one can see that the effect of large S in the exponent is suppressed by a a small factor (/l/lSl) 1 / 2 . There are 
two ways to achieve this suppression and to increase A. The first is to use magnetic clusters with very strong easy 
plane anisotropy and relatively weak easy axis anisotropy in that plane. Particles of f£h and Dy may satisfy this 
condition. The second way is to place the particle above the surface of a superconductoill3. In that case the magnetic 
dipole interaction of S with its mirror image inside the superconductor effectively reduces the uniaxial anisotropy 
A. The particle and the superconductor should be selected such that the magnetic field induced by the particle, 
H ~ AitubS, does not exceed the first critical field of the superconductor, H c \. Manipulating the distance between 
the particle and the superconductor, one can achieve the condition A « \B\. This can be done by, e.g., controlling 
the distance with a ferroelectric buffer in the electric field. Note that in the absence of the external magnetic field the 
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odd symmetry of S with respect to time reversal preserves the coherence of such a setup independently of the shape 
of the superconducting surface and electric fields in the problem. 

The tunneling of S also can be induced by its hyperfine interaction with nuclear spins, V = BS1, where I = J^Ij 
is the total nuclear spin of the clusterj-obtained by summing over spins of individual nuclei. This problem is rather 
involved. It has been studied in RefJl3 and is relevant to tunneling in Mn-12. The total Hamiltonian conserves 
the magnitude of I and the Z-projection of S + I. In the millikclvin range nuclear spins must order, developing 
Imax = y~)\Ij\- It is easy to see that the problem is the one of quantum coherence only if I m ax = S. In that case 
the classical ground states correspond to S and I of equal length looking opposite to each other along the Z-axis, 
\S > | — Imax > and I — S > \I ma x >■ Tunneling removes the degeneracy of the ground state. The corresponding 
splitting can be obtained by the perturbation theory for B << A: 



A = 8(A + B)S 2 
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In Mn-12 S = 10 while Imax — 30, so that the coherence of the above type is impossible. It is not out of the question, 
however, that chemists will produce a molecular cluster with S = / in the future. 



B. Antiferromagnetic clusters 

Tunneling in antiferromagnetic clusters0 turns out to be much stronger than in ferromagnetic clusters, making them 
promising candidates for quantum coherence. Consider an anisotropic antiferromagnetic cluster with two compensated 
sublattices of spin Si and S2, described by the Hamiltonian 

H = -A(S 2 lz + S 2 2z ) + BS r S 2 (9) 

with positive A and B satisfying A « B. 

Let us show that this rapdel can be mapped onto the model with strong transverse anisotropy. The Lagrangian 
corresponding to Eq.(^J) iscl 



L = S((f>i cos 6»i + 4> 2 cos6> 2 ) - S(<j>i + <j) 2 ) + AS 2 {cos 2 9 1 + cos 2 9 2 ) - BS 2 cos8icos8 2 - BS 2 sm9 1 sm9 2 cos(4>i 
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where <p±, 9±, <fi 2 , 9 2 are spherical coordinates of vectors Si and S2 of fixed length S. The path integral over these 
angles is dominated by cos6>i = — cos#2 = cos#, sin^i = sin6> 2 = sin6>. Introducing <fi = <\>\ — 2 , one obtains, up to a 
phase term, the effective Lagrangian 

L eff = S<j>cos6 + (2A + B)S 2 cos 2 6 - BS 2 sm 2 9cos<f> . (11) 

With the notations 0/2 = $ and 25 = <t, it can be transformed into 

A B 

L e ff = a$cos9 + —a 2 cos 2 9 - —<j 2 sm 2 9cos 2 <& , (12) 



which is equivalent to the Hamiltonian 



H = —o 3 , + f<7* • (13) 



One can then use the known result, Eq.(7), to obtain the tunneling splitting, 
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Here B is the exchange constant which is typically 10 4 — 10 6 the anisotropy constant A. Consequently, antiferromag- 
netic particles consisting of a few thousand magnetic atoms can exhibit a significant tunneling rate between the |ti > 
and ||t > states. 

One should notice, that, in order to manipulate these states by the magnetic field, some magnetic non-compensation 
of the sublattices is needed. In this case tunneling is possible only due to the presence in the Hamiltonian of the 
transverse field, or the transverse anisotropy, e.g., b(S 2 x + S 2 X ). Let the non-compensated spin be s. It has been 
demonstrated!^ that the antiferromagnetic tunneling with A oc exp [-A-SiA/B) 1 / 2 ] holds upto s^^/Bf^S « S. 
At greater s it switches to the ferromagnetic tunneling with Aoc exp[— 2s{A/b) 1 / 2 ]. Strongly non-compensated ferri- 
magnetic clusteijSj-Jike, e.g., Mn-12, are always in the ferromagnetic tunneling regime, while weakly non-compensated 
ferritin particlesQ'B can be in the antiferromagnetic tunneling regime. 
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III. DECOHERENCE 



A magnetic cluster of the type described above will always be imbedded in a non-magnetic solid dissipative environ- 
ment. The potential for low decoherence arises from a number of reasons. The first of them is that strong electrostatic 
interactions are involved, through exchange couplings, only in the formation of the single spin S of the cluster. All 
other interactions of S have relativistic smallness of order (v/c) 2 to a some power. Due to this fact the ferromagnetic 
resonance in some materials has a quality factor of one million. The second reason is selection rules for spins due to 
the time reversal symmetry discussed below. 

In quantum computation one is interested in creating an arbitrary superposition of the | f > and | J_> states, 

\*>=C l \1>+C 2 \l> . (15) 
Using equations (2) and (3) this state can be re- written in terms of |0 > and |1 >: 

|* >= C[\0 > +C' 2 \1 > , (16) 

where C[ = {C\ + C 2 )/V^ and C 2 = {C\ — C 2 )/V2- It is clear, therefore, that the spontaneous decay of the excited 
state |1 > into the ground state |0 >, accompanied by the emission of the energy quantum A, should be a major 
concern for preserving quantum coherence. 

Note that, in principle, there may be decohering processes involving other excited levels of the spin Hamiltonian. 
For large spin, due to a small tunneling rate, such levels are separated from |0 > and |1 > by the energy gap (say A) 
that is large compared to A. An example of such a process would be an Orbach two-phonon process corresponding to 
the transition |1 > — > \A > caused by the absorption of a phonon, followed by the spontaneous decay \A > — > |0 > 
with the emission of a phonon. Also, there may be processes involving the excited states of the environment. An 
example of such a process would be a two-phonon Raman process that corresponds to the emission and absorption of 
two real phonons satisfying h(uji — lu 2 ) = A. All such processes are strongly temperature-dependent. Theix-strength 
is measured by exp(— A/ksT) or by some high power of T/G where @~1Q 2 K is the Debye temperaturaLj. In the 
millikelvin range the rate of such processes is negligible. Consequently, they are of little concern for the decoherence. 
On the contrary, processes of the spontaneous decay |1 > — > |0 > can exist even at T=0. One should therefore 
concentrate on such processes. 

As follows from the previous section, cases of interest for quantum coherence are described by Hamiltonians which 
contain even powers of spin operators. Due to the time reversal symmetry, this will always be the case in the absence 
of the external magnetic field. Consequently, the states |0 > and |1 > most generally can be written as 

s 

|0> = J2 a ™\ m > ( 17 ) 
m=-S 
S 

11 > = 53 /3 ™' m > ' 

m=~S 

where a m — ot- m and (3 m = — /3_ m . These equations are exact, as compared to the approximate equations (2) and 
(3). They simply reflect the fact that \m > is a complete set of vectors in the Gilbert space of the spin Hamiltonian 
and that |0 > and |1 > have different symmetry with respect to time reversal. 

Let K be the antilinear antiunitary operator of time reversal. The spin operator is odd with respect to time reversal, 
KSK^ = -S. On the contrary, the spin Hamiltonian that only contains even powers of components of S is even with 
respect to time reversal, KHK^ = H . Consider now a decohering operator D. D can be due to the interaction of S 
with phonons, electromagnetic fields, nuclear spins, etc. These interactions have different symmetry with respect to 
time reversal. Let D Q and D e be time-odd and time-even operators respectively. Since the state \0 > is even with 
respect to time reversal and 1 1 > is odd, the following general statement is true 

<0|O e |l>=0 . (18) 

The spin-phonon interaction is of the form 

H sp — aikirnSiSk^—^- + h.c. , (19) 

where u is the lattice displacement and aikim is a tensor reflecting the symmetry of the lattice. In terms of the 
operators of creation and annihilation of phonons 
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T, T:^i/2 ^- a ix) > (20) 



(2MiV)V2 (wfcA )i/a 

where M is the unit cell mass, 2V is the number of cells in the lattice, ej-.A is the phonon polarization vector, A = t, t, I, 
and u>k\ — v\k is the phonon frequency (v\ being the speed of sound). The spin-phonon interaction given by Eq.(19) 
describes transitions between different spin states accompanied by the emission and absorption of phonons. For 
instance, in the absence of tunneling the excited state \S — 1 > of the Hamiltonian H = —AS 2 , corresponding to tha 
ferromagnetic resonance, can relax to the ground state IS* > at T = by spontaneously emitting a phonon at a ratecJ 

3 



A 2 <> 2 A 
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where p is the mass density of the lattice and C is a constant of order unity. For, e.g., A~ IK, S*~ 10 3 , (A/7i) ~ 10 9 s _1 , 
p~ lg/cm 3 , and v~ 10 5 cm/s, Eq.(21) gives 7~ 10 3 s _1 . 

We should now notice that KH sp K^ — H sp , that is, the spin-phonon operator is even with respect to time reversal. 
Consequently, < 0|if sp |l >= 0. This also can be checked by the direct calculation of matrix elements of H sp using 
expressions (17). Thus, |1 > cannot spontaneously decay into |0 > with an emission of a phonon and we conclude that 
in a millikelvin range the spontaneous emission of phonons cannot decohere the coherent superposition of the | |> and 
| {> spin states. This is a strong statement which requires some clarification. Indeed, the spin-phonon interaction 
originates from the spin-orbit coupling of the form iJ so ocL-S where L is the orbital momentum. Its density can 
be presented as petkirkUi- This operator seems to have non-zero matrix elements between |0 > and |1 > given by 
equations (17). However, this is only because the formulation of the spin tunneling problem presented above does not 
account for the conservation of the total angular momentum. In fact, the coherence is possible only if L + S = OcJ. 
Consequently, as S tunnels between | "f> and | |>, so does the mechanical angular momentum of the body. This 
is an analogue of the Mossbauer effect for spin tunneling. The energy associated with the mechanical rotation is 
h 2 S 2 /2Ii„ where Ii n is the moment of inertia of the solid matrix containing the magnetic cluster. This energy should 
not exceed A, otherwise it would be energetically favorable for S to localize in the | |> or | i> state. Since Ii n scales 
as the fifth power of the size of the matrix, it is easy to see that such localization may occur in a free particle of 
size less then 10 nm, while for bigger systems the conservation of the total angular momentum is rather formal than 
practical question. With the account of the momentum conservation, the ground state and the first excited state are 
> | — L > ± | — S > \L >), where the absolute values of S and L are equal. Now, again, the |0 > state is even 

with respect to time reversal while the |1 > state is odd. The spin-orbit operator is even, KH so K^ = H so , because L 
is odd, KhK^ = -L. Thus, as before, < 0|if so |l >= 0. 

As has been shown in the previous section (see also Refs. 21,22), nuclear spins always destroy the coherence at 
H = unless tunneling is induced by the hypcrfine interaction with I = S. Except for that exotic possibility, nuclear 
spins must be always eliminated from magnetic qubits by the isotopic purification. Similarly, the presence of free 
non-superconducting electrons in the sample will decohere tunneling through the spin scattering of electrons, D e cx s-S. 
Although, this operator is time-even, free electrons incidentally passing through the magnetic cluster, perturb |0 > 
and |1 >, breaking their properties with respect to time reversal. Thus, strongly insulating materials should be chosen 
for magnetic qubits. The effect of incidental phonons due to, e.g., relaxation of elastic stresses in the matrix (the 1// 
noise), should be similar to the effect of incidental electrons in perturbing |W >. Thus, the perfection of the lattice 
should be given a serious thought when manufacturing magnetic or any other qubits. One should then worry about 
the decohering effect of spin interactions which are odd with respect to time reversal. These are Zeeman terms due 
to magnetic fields. For example, D D = —gp,BS z H z (t) has a non-zero matrix element between |0 > and |1 >. The 
effort should be made, therefore, to shield the magnetic qubit from the magnetic fields during the process of quantum 
computation. This can be done by placing the magnetic cluster inside a nanoscopic superconducting ring. Such a 
ring may be used to control and measure the states of the qubit. Connecting such rings by superconducting lines may 
be the way to make a miltiqubit system. 
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